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Abstract 

We analyze the possibility that the lighter stop t\ could be the next-to-lightest super- 
symmetric particle (NLSP) in models where the gravitino is the lightest supersymmetric 
particle (LSP). We do not find any possibility for a stop NLSP in the constrained MSSM 
with universal input soft supersymmetry-breaking masses at the GUT scale (CMSSM), but 
do find small allowed regions in models with non-universal Higgs masses (NUHM). We dis- 
cuss the cosmological evolution of stop hadrons. Most t\qq 'sbaryons' and the corresponding 
'antisbaryons' annihilate with conventional antibaryons and baryons into tiq 'mesinos' and 
the corresponding 'antimesinos', respectively, shortly after the quark-hadron transition in 
the early Universe, and most mesinos and antimesinos subsequently annihilate. As a result, 
insufficient metastable charged stop hadrons survive to alter Big Bang nucleosynthesis. 
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1 Introduction 



In many supersymmetric models there is a multiplicatively-conserved quantum number, R 
parity, that guarantees the stability of the lightest supersymmetric particle (LSP). In order 
to avoid the LSP binding to ordinary matter, it is usually assumed to have neither strong 
nor electric charge [1]. Candidates for the LSP in the minimal supersymmetric extension of 
the Standard Model (MSSM) with gravity include sneutrinos, the lightest neutralino x an d 
the gravitino G. Light sneutrinos were excluded by searches for invisible Z decays at LEP, 
and heavier stable sneutrinos would have been found in direct searches for the scattering of 
astrophysical dark matter particles on ordinary matter [2] . Thus, most attention has focused 
on the neutralino and the gravitino. Overlooked to some extent, a gravitino LSP is in fact 
quite generic even in models based on minimal supergravity (mSUGRA) [3]. 

In the case of a gravitino LSP [4-10], the next-to-lightest supersymmetric particle (NLSP) 
has a long lifetime, decaying with gravitational-strength interactions if supersymmetry break- 
ing is mediated by supergravity. The question of the identity of the NLSP then becomes 
important. One generic possibility is that the NLSP is the lightest neutralino, in which case 
the long-lived x would probably decay unseen, mainly via x ~^ C + 7, without being stopped 
beforehand [4,5]. Another generic possibility is the lightest charged slepton, probably the 
lighter stau fi in the MSSM with universal scalar soft supersymmetry-breaking masses (the 
CMSSM) [4,6,7,9]. This leads to scenarios with a metastable charged sparticle that would 
have dramatic signatures at colliders [6,11-14] and could affect drastically the cosmological 
abundances of light elements. 

Electromagnetic showers from the decay products of metastable particles can alter the 
abundances of light elements by photo-dissociation and subsequent secondary reactions [15- 
17]. Moreover, hadronic showers can alter the amounts of baryons involved in the Big-Bang 
nucleosynthesis (BBN) processes if the lifetime ^ 10 6 s [18]. However, a more significant 
effect can occur in the case of a negatively- charged particle, which can form an electromag- 
netic bound state with a nucleus, and influence the BBN processes by lowering the Coulomb 
barrier for nuclear fusion [19,20] (a catalytic effect). This has been studied within the GDM 
in the case of a stau NLSP in some CMSSM and mSUGRA scenarios [21]. 

However, there are also other possible candidates for the NLSP, such as some sneutrino [7, 
22] or squark species. Among the different squark species, a generic candidate for the lightest 
is the lighter stop t\ [23] , which would have interesting implications for cosmology [24, 25] , 
although there are other possibilities. In this paper we study the feasibility of scenarios 
with a gravitino LSP and a stop NLSP. Thus, we search for regions of the MSSM parameter 
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space where the t\ is lighter than the supersymmetric partners of all the other Standard 
Model particles, including the neutralino x- A previous study showed that this is possible 
for large values of the soft trilinear supersymmetry-breaking parameter Aq, and the regions 
of the CMSSM parameter space where this happens have been delineated [25]. In traditional 
scenarios with conserved R parity and a heavy gravitino, these regions would have been 
discarded because they have a charged and coloured stable particle. However, if the gravitino 
G is the LSP and therefore constitutes the dark matter, one should explore whether some 
parts of these regions might survive. 

There are several experimental and cosmological constraints on such a stop NLSP sce- 
nario that must be taken into account. As discussed in more detail below, the lifetime of 
the t\ may be (very) long, in which case the relevant collider limits are those on (appar- 
ently) stable charged particles. We interpret the limits available from the Tevatron collider 
as implying that > 220 GeV [26] □. We find no regions of the CMSSM parameter 
space compatible with this and other experimental constraints on t\ NLSP scenarios. How- 
ever, when we relax the CMSSM universality assumptions by considering non-universal soft 
supersymmetry-breaking masses for the Higgs fields (NUHM models), we do find limited 
regions of parameter space with a t% NLSP. Typical allowed values of the NUHM parameters 
are m 1/2 ~ 600 GeV, m ~ 500 GeV, A ~ 2100 GeV, n ~ 750 GeV, m A ~ 1400 GeV and 
tan/3 ~ 10. 

We then consider the cosmological constraints on such cases. As we show, the density of 
t\ sparticles and antisparticles after cosmological freeze-out at a temperature of several GeV 
is strongly suppressed by strong couplings in the annihilation processes. Subsequently, at the 
quark-hadron transition these stops would have combined with quarks into t\qq 'sbaryons' 
and t\q 'mesinos' and the corresponding antiparticles. The late decays of these stop hadrons 
could have affected the light-element abundances obtained from Big-Bang nucleosynthesis, 
and negatively- charged antisbaryons and antimesinos could have had dramatic bound-state 
effects. However, we argue that the great majority of the stop antisbaryons would have 
annihilated with conventional baryons to make stop antimesinos, and that most mesinos and 
antimesinos would subsequently have annihilated [27]. Any negatively-charged antimesinos 
would have decayed (relatively) rapidly into neutral mesinos. These would have been (al- 
most) the only metastable t\ relic particles, and would be relatively innocuous, despite their 
long lifetimes, because they would not have important bound-state effects. Because of the 
low density of t\ after freeze-out following coannihilation and the subsequent cosmological 
evolution, this limited region of stop NLSP scenarios within the NUHM framework seems 

1 The LHC will probably be sensitive to a metastable t\ that is an order of magnitude heavier. 
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to be viable. We conclude our paper with a brief discussion how such a scenario could be 
probed experimentally. 



2 Stop Properties 

2.1 Stop Masses, Mixing and Couplings 

We start by giving some important formulae and making some crucial definitions. The (2x2) 
stop mass matrix may be written as: 



Ml 



( Ml L M\ R ^ 

Mil M RR 



(1) 



where the entries take the forms: 

M 2 LL =M? L +m 2 t +l cos 2(3 (4m 2 w - m|) , 

M 2 RR = M 2 r + m 2 + | cos 2(3 sin 2 6 W m% , (2) 
m lr = -mt{A t + n cot (3) = -m t X t . 

The corresponding mass eigenvalues are given by: 

ml=m 2 + ^(Ml+Ml) + ^m 2 z cos2(3-^, (3) 

and 

ml = m 2 + \(Ml+ M?J + X -m\ cos 2(3 + -, (4) 

where A 2 = (m~ l — M 2 ^ + | cos 2(3{Sm 2 v — hm 2 z )\ + Am 2 \A t + fi cot (3\ 2 . The mixing angle 
Ql between the weak basis and the mass eigenstates (ii,^), is given by tan# t ~ = 

(m~ — M 2 L )/\M 2 R \. It is clear that obtaining a very light stop requires a very large value 
for the trilinear soft supersymmetry-breaking parameter [25]. 

The interactions of the left and right antistops t* L and t* R with the gravitino field ^ and 
the top quark are given by [28] B 

c = --^[^,rrp R tdj R + *^p L td v t* L \, (5) 

where the reduced Planck mass is given by M = M p i/^/8tt, with M p i = 1.2 x 10 19 GeV. The 
interaction lagrangian for t\ t 2 is then: 

C = --^[^7V(sin^ + cos6 { P L )tdJl + ^ 7 V(cos%P R - sin 6 { P L ) t d u t* 2 }. (6) 



2 Note, however, that there is a typographical error in Eq.(4.31) of Ref. [28]: {^—m 3 / 2 ) should be (jji+m^^)- 
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The corresponding Feynman rule for the vertex is: 

H (p) & t - - -i^ rfi (sin % P R + cos 9~ t P L ). (7) 
Similarly, the Feynman rule for the chargino-gravitino-iy vertex is: 

XT&W-ik) - -^7Y(4,P fl + A m P L ) (8) 

Here A Li = U* 2 cos /3, and A Ri = V* 2 sin /3, where K and Z7 are the matrices that diagonalize 
the chargino mass matrix. 

2.2 Stop Decay Modes and Lifetime 

There are several possible scenarios for stop decay, depending on the mass difference between 
the stop NLSP and the gravitino LSP Am = raj ~ m ci anticipating that a stop NLSP must 
have m^ > m t from the direct search bound. 

1. Case 1: Am > m t , i.e. small mg ^ m^ — m t . In this case, the stop can decay directly 
into a top quark and a gravitino, and the rate for this dominant decay is 

T = — k — t- [4 (m~ — m~ — mV\ + 20 sin 9t cos 9jm t mp\ 

192tt Mlm%m 3 r 1 v * x u * * Gj 

+ TrV~ — ) 2 — 4m^m~] [( m ^ + m t — m g) 2 — 4m? i m^] ^ 2 . (9) 



x 



This decay rate is similar to that for stau decay into tau plus gravitino [4] , but in this 
case m t cannot be neglected. Previous results are reproduced in the limits m t — > and 
0~t - o. 

We show in Fig. [TJsome typical numerical results for m^ = 200, 300, 400 and 500 GeV 
(from top to bottom), mg < mj — m t and m t = 171 A GeV [29], for both zero stop 
mixing (red solid line) and maximal mixing (blue dashed line) We see that the stop 
lifetime is relatively insensitive to the stop mixing angle but depends sensitively 
on the sparticle masses, and ranges between 10 3 and 10 9 s. Clearly, this is extremely 
long compared with the QCD hadronization time-scale, so that the stop NLSP (unlike 
the t quark) forms metastable hadrons, whose spectroscopy and phenomenology we 
consider below. Moreover, this lifetime range is also very long on the typical time- 
scales of collider experiments, which must therefore consider how to search for these 
stop hadrons. Indeed, the stop lifetime fits into the range where the cosmological 
effects considered in [16,18,19,21] become important. 



3 Typical values in the allowed stop NLSP region in the NUHM are 6i ~ 1.3. 
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m 5 (GeV) 

Figure 1: The stop lifetime as a function of for = 200,300,400 and 500 GeV 
(top to bottom), shown for the case of the two-body decay t\ — > Gt, the dominant mode for 
< —m t , assuming zero stop mixing (red solid line) and maximal mixing (blue dashed 
line). 

2. Case 2: my/ + mj < Am < m t . In this case, the dominant decays are into the 
three-body final state t\ — > G + W + b. We identify three tree-level decay diagrams, 
proceeding via t, b and chargino exchange. The amplitudes are 

M t = C t P t (q 1 )^^[A i + B r y 5 ](4 1 + rn t )Ye P (k)P L v(p 2 ), (10) 
M~ hi = C^q^.q^PL + kPnjp^vipi), (11) 
M x + = C , x+ P Xi+ (g 3 )^7%(^)7' 1 ^ + ^75](^ + m x )[^ + P,75]^(P2), (12) 

where C t = g 2 /2M, Cg. = 2g 2 Ki/M, C x + = mw/M, and p is the initial stop four- 
momentum. We define q± = p — pi, q 2 = p — k and q 3 = p — p 2 , with pi, k,p 2 denoting 
the outgoing four-momenta of the gravitino, W boson, and b quark respectively, and 
e p (k) denotes the W polarization vector. Expressions for Af, Bf, a i; b i: Vi, A { , Si and 
Pi are presented in the Appendix. 

Squaring and summing over final polarizations we obtain: 

\M\ 2 = \M t \ 2 + \M- b \ 2 + \M x+ \ 2 + 2Re[M\M~ b + M\M X+ + M\M x+ }. (13) 
where the sums over sbottom and chargino indices are implicit. The individual squared 
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amplitudes can be written as: 

\M^\ 2 = Cl\PM\ 2 W Ma , (14) 

where ipa = (t, bj, xt)i and the functions W^ a ^ a are functions of the scalar products of 
the momenta p,Pi,P2, k. The interference terms may be written as follows: 

M\M^ b = c;c^ b p; a (q a )PMw Mb , (is) 

where the functions W$ a $ b can be written also in terms of the invariants. The functions 
Pipaila) are propagator factors, e.g., for the top quark ip a = t, and we have 

PMi) = - — ( 16 ) 

q{ — m t + le 

There are similar expressions for the sbottom and chargino contributions, P b {q2) and 
P x +(q 3 ) respectively. 

Detailed formulae for the functions W^, a ^ a and W^ a1 p b are given in the Appendix. Using 
these, we calculate the decay width: 

dV - ml '.M| 2 , (17) 



dx dy 256 tt 3 

where the integration limits are (in the limit when we neglect the bottom quark mass): 
2/ic < x < 1 + fiG — and ?/_ < y < y + , where Ui = mf/mj and: 

^T 

^+^G + ^W-X W . , , 1/2l , , 

2(1 + hg-x) 

Integrating this equation numerically, we find that in significant regions of parameter 
space the light stop has a lifetime of order 10 9 — 10 14 s or more. Typical results are 
shown in Fig. As expected, the typical lifetimes are orders of magnitude longer than 
those in Fig. [TJ They are also sensitive mainly to the sparticle masses, and relatively 
insensitive to the stop mixing angle, as can be seen by comparing the red solid and 
blue dashed curves, as well as being insensitive to the sbottom mixing angle 9~ bl which 
is assumed here to vanish 0. 

3. Case 3: m b + Aq C d < Am < m w + m b . In this case, the real W of the previous case 
must become virtual, and the dominant decays are four-body: t\ —>■ G + b + qq or iv. 
The decay rate for this case is further suppressed compared to Case 2, and we estimate 



4 Typical values in the allowed NUHM models discussed later are Or ~ 0.04. 
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Figure 2: The stop lifetime as a function of tUq for m~ tx = 200, 300, 400 and 500 GeV, shown 

for the case of three-body decays t\ — > GWb, the dominant mode for mw + mb < Am < m t , 
assuming zero stop mixing (red solid line) or maximal mixing (blue dashed line). 

that the stop lifetime would in this case exceed 10 12 s. Thus, the stop might even decay 
after the release of the CMB, in which case there would be important constraints from 
the absence of distortions in CMB spectrum. We have not explored this issue, in view 
of the relatively small region of parameter space concerned. 

It is apparent from the above discussion that not only does the stop live long enough 
to hadronize and pass through collider detectors, but it may also live long enough to wreak 
cosmological havoc. We discuss each of these aspects in the following Sections. 

3 Spectroscopy of Stop Hadrons and their Decays 

The metastable stop would hadronize to produce both t\qq 'sbaryons' and t\q 'mesinos' and 
their antiparticles, many of whose aspects are discussed in [30]. On general QCD principles 
and by analogy with the spectroscopy of charmed hadrons, one expects the T° = t\u, T + = 
t\d and T s =t\s mesinos to be the lightest stop hadrons. As was pointed out in [30], one can 
expect the T° mesino and its antiparticle to be strongly mixed. Since m s > > m u , and 
since the T + and T s mesinos, being charged, would acquire additional electromagnetic mass 
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corrections, we expect mf a — mf + and mf+ — rrif to be similar to the measured values of the 
D s — D + and D + — D° mass differences, namely ~ 99 MeV and ~ 4.8 MeV, respectively [31]. 

Correspondingly, we would expect the T s mesino to decay weakly into T°ez/ with a lifetime 
similar to that of the muon, namely ~ 2 x 1CT 6 s, and the T + mesino to decay weakly into 
T°ev with a lifetime 

/ m n -m p V 1 o /,m 

Tf+ ~ T n x — ~ 1.2 s. (19) 

\rrif + - m^ / 

These mesino lifetimes are also such that they would pass through a typical collider detector 
before decaying. In the early Universe, the T s would have decayed very quickly after being 
formed at the quark-hadron transition, whereas the T + , if they survive, would have decayed 
near the beginning of BBN, and so would not have affected its end results. Moreover, since 
the stable T° mesino would be neutral, it could not have catalyzed light-element nucle- 
osynthesis by bound-state effects. The only potential cosmological danger from the mesinos 
would be the supersymmetric decays of the T° into the G and conventional particles, as 
discussed in [21], for example. 

Turning now to the t\ sbaryons, again by analogy with the charmed baryons, we would ex- 
pect the lightest state to be the A~ = tiud, with the other sbaryons = ti(uu, ud, dd), 
S~'° = tis(u,d) being heavier by amounts ~ Aqcd^s — n^d,u, respectively. Just like the 
T s mesino discussed above, these heavier sbaryons would have decayed innocuously before 
BBN. For example, if the — mass difference were similar to the corresponding mass 
differences among charmed and bottom baryons, namely ~ 170 [31] to 190 MeV [32], the 
would decay very rapidly via the strong interactions. If the Ef — Af mass difference were sim- 
ilar to the corresponding mass difference among charmed baryons, namely ~ 180 MeV, the 
S~ = tisd state would decay semileptonically with a lifetime < 10~ 6 s into the A^, whereas 
the S~ = tisu state would decay semileptonically into the S~ + = t±uu with a longer lifetime 
~ 10~ 2 s (because of the much smaller phase space ~ 15 MeV for the decay). However, 
this decay would also occur by the beginning of BBN. Thus, the cosmological dangers could 
potentially arise only from the supersymmetric decays of the A~ and its antiparticle, the A~, 
and dangerous catalysis effects could only be due to bound states of the A~. 

We note in passing some similarities with and differences from the case of a metastable 
charge -1/3 squark, such as the lighter sbottom, b±. In this case, we would expect the 
lightest sbaryon to be the A~ = b\ud. Its decays might cause cosmological problems, but it 
could not cause dangerous bound-state effects. The nature of the lightest sbottom mesino 
is not so clear. The fact that the d quark is heavier than the u quark would tend to 
make the B° mesino heavier, but the electromagnetic corrections would add to the mass 
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of the B + . Experimentally, the situation with B mesons is ambiguous, m B o -m B + = 
0.33 ±0.28 MeV. However, the likelihood is probably that the B + would be lighter, in which 
case its antiparticle, the B~ , would generate bound-state effects, unlike the T°. 

4 Collider Lower Limit on the Stop Mass 

As we have seen, the stop would have a very long lifetime in GDM scenarios of the type 
considered here, and would hadronize before passing through a typical collider detector. 
The relative production rates of mesons and baryons containing heavy quarks are not es- 
tablished, and neither are the relative production rates of heavy-quark mesons containing 
strange quarks. We assume for simplicity that half of the produced stops hadronize into 
charged mesinos or sbaryons, and half into neutral stop hadrons. These would be produced 
embedded within hadronic jets, but conventional QCD fragmentation ideas suggest that the 
stop hadrons would carry essentially all the energies in these jets, with energy fractions 
Zf~l- K QCD /m~ tl . 

The typical energy loss as the stop hadron passes through a detector tracking system 
would be very small. There would also be nuclear interactions, particularly in calorime- 
ters. In addition to the familiar charge-exchange reactions, these would also include baryon- 
exchange reactions, whereby a stop mesino striking a nucleus would convert into a stop 
sbaryon: T + (p, n) — > (A^, S^) + nn, whereas the corresponding sbaryon-to-mesino con- 
version would be essentially forbidden. It has been pointed out that the baryon-exchange 
process would be almost 100 % efficient in converting heavy mesinos to sbaryons when they 
traverse material with a thickness of 1 m or more [30] (see also [33]). We therefore assume 
for simplicity that all of the stop hadrons emerging from calorimeters into muon detectors 
are sbaryons, that half of them are singly-charged, and that this charge is independent of the 
charge of the stop hadron at production. This would imply that just a quarter of the pro- 
duced stop hadrons would be singly- charged both at production and in the muon detectors. 
Thus, only about 1/16 of the produced stop-antistop pairs would yield a robust signature of 
a pair of oppositely-charged massive metastable particles. 

We use here the limits set by direct searches for the pair-production of massive (meta)stable 
charged particles at hadron colliders to set a lower bound on the stop mass. Nunnemann [34] 
gives an upper limit from DO of about 0.1 pb from a search for the pair-production of massive 
oppositely-charged particles, and a similar upper limit for the pair-production of stops has 
been presented by CDF [26], which is used to set a lower limit of 220 GeV for m^. Gallo [35] 
gives a DO upper limit on the production of neutral gluino hadrons of about 0.5 pb. Again 



9 



assuming that about a quarter of the stop hadrons are produced neutral and also appear 
neutral in the outer detectors, this limit gives a somewhat weaker limit on vn- tl . Therefore, 
we assume > 220 GeV [26]. 

5 Stop NLSP in the CMSSM 

We now discuss the prospects for finding a stop NLSP in the CMSSM, i.e., the simplest 
variant of the MSSM, in which the soft supersymmetry-breaking masses are universal at the 
GUT scale. Thus, we have as free parameters mo, the universal soft scalar mass at the GUT 
scale, mi/2, the universal gaugino mass at the GUT scale, A , the universal trilinear soft 
supersymmetry-breaking parameter at the GUT scale, tan/3, the ratio of the two MSSM 
Higgs vevs, and the sign of // (where /i is the Higgs mixing parameter). In addition, unless 
we make additional assumptions, as may be motivated by supergravity models, we must 
consider the gravitino mass mg as an extra free parameter, which is chosen so that the 
gravitino is the lightest supersymmetric particle (LSP). 

We search within the CMSSM for a set of parameters where not only is the stop the NLSP, 
but also all the known experimental and phenomenological constraints on supersymmetry 
are satisfied, including the b — > 57 decay rate [36,37], and the LEP lower bounds on the 

B 



masses of the chargino and the Higgs boson [38] □. In view of the ambiguity in the value 
of the hadronic contribution to the Standard Model value of the muon anomalous magnetic 
moment, — 2, we omit this observable for the moment. The B s — > fi + fi~ constraint is 
significant only for large tan /3, whereas, as we see below, the regions that are relevant to our 
search have relatively small tan/3 ~ 10. 

We choose the sign of /i to bepositive, as our search indicates that negative /i has less 
chance of yielding a stop NLSP and assume m t = 171.4 GeV [29] and mb(mb) MS = 
4.25 GeV. Since m x ~ 0.43m!/ 2 , in order to obtain m x > ntf > 220 GeV, we must set 
m i/2 ^ 520 GeV in the CMSSM (this is relatively independent of tan/?). However, the 
LL and RR components of the m~ mass matrix receive contributions of about 6m^ 2 and 
4m^ 2 respectively, forcing one to consider large off-diagonal elements. These are of the 
form m t X t = —m t (A t + [ij tan/3), as seen in (2) Low stop mass eigenvalues therefore 
require a combination of high values of A§ and relatively low values of tan/3. In fact, as 
we show explicitly below, only intermediate values of tan/3 have any chance of realizing 



5 We use the public Fortran code FeynHiggs [39] to calculate m/,. 
6 Negative \i may also be disfavoured by b — > 57 and g M — 2. 

7 Note that A t signifies the value of the trilinear term at the scale , which differs from its value Aq at 
the GUT input scale. 
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low stop masses without upsetting the remaining phenomenological constraints. At high 
tan/? the b — > s'j constraint is not satisfied, and at low tan/3 the LEP Higgs mass bound 
cannot be satisfied. Our only option therefore is large Aq. Our search in the CMSSM is 
further complicated by the dependence of the Higgs mass on X t . For relatively small values 
of \X t \/mi/2, the Higgs mass increases with increasing |Jf t |. However, for \X t \/mi/ 2 > 2, 
the Higgs mass begins to decrease rapidly [40]. In order to obtain a light stop, we need 
\X t \/mi/2 ~ 6mi/ 2 / m t, which is too large to have any chance of satisfying the Higgs mass 
constraint. 

In order to obtain a more comprehensive picture, we show in Fig. [3] some contour plots in 
(tan/?, Aq) planes for some fixed values of mo and m\/2- In panel (a) we fix itlq = 500 GeV 
and mi/2 = 400 GeV (the gravitino mass is irrelevant here). The b — > 57 constraint, 
which excludes large tan/?, is shown by the green shaded region. In the allowed region, 
either the stop or the neutralino is the lightest sparticle in the spectrum. Above the thick 
purple solid line, the stop is lighter, while below it the neutralino is lighter. We also plot the 
mi = 220 GeV contour, represented by the orange solid line: regions above this line have 
m~ tl < 220 GeV, and therefore are excluded. The Higgs mass constraint is represented by two 
red lines, the dashed line is based on a likelihood analysis, and the dot-dashed line on the face 
value of the Higgs mass limit, namely = 114.4 GeV (deprecated). The constraints are 
satisfied only below the lines. We use the constraint determined by the likelihood function 
in our analysis. We see in panel (a) that there is an overlap between the stop NLSP region 
and the region allowed by the Higgs likelihood constraint. However, the stop mass is around 
150-160 GeV in this region, which is therefore excluded by the lower bound on the stop 
mass. We also plot in Fig. [3](a), the contour where the NLSP would have a relic density of 
Q^h 2 = 4 x 10 -4 , if it did not decay. This is shown by the thin green line that lies below 
the neutralino-stop contour. The small value is a consequence of the strong stop-antistop 
annihilation cross section. 

When we increase mi/2, for example to 450 GeV as in panel (b), both m x and raj 
are raised and neutralino-stop degeneracy is reached at higher Aq. We see that the m^ = 
220 GeV line is closer to the neutralino-stop degeneracy line, and if we kept increasing 
we would be able to find points where the stop is lighter than the neutralino and has mass 
larger than 220 GeV. However, although the Higgs mass constraint also moves to higher Aq, 
it moves slower than the previous two lines. As a result, there is no overlap region where 
there is a stop NLSP and the Higgs mass bound is satisfied. 

The first two panels in Fig [3] already suggest that there are no allowed regions with a 
stop NLSP and gravitino LSP in the CMSSM. Generalizing this observation, we first note 
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Figure 3: The (tan/3,A ) p/ane zn Sfce CMSSM, for (m 1/2 ,m ) = faj (400,500) GeV, (b) 
(450,500) GeV, (c)j400,450) GeV, and (d) (400,550) GeV respectively. We use m t = 
171.4 GeV, rrib(mb} MS = 4.25 GeV, and /i > 0. The neutralino stop degeneracy contour is 
plotted as the thick solid purple line: above this line stop is lighter and the NLSP, assuming 
a light gravitino. The solid orange line is the contour for = 220 GeV: above this line, 
the stop is too light. Large tan/3 is excluded by the b —+ 57 constraint (green shaded region). 
The Higgs likelihood exclusion line (preferred) is drawn as a dashed red line, while the face 
value rrih = 114.4 GeV (deprecated) is the dot-dashed line. Also plotted as the thin green line 
is the stop relic density Qjh 2 = 4 x 1CT 4 . 
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that, due to the nature of the RGEs in the MSSM, varying mo would have less effect than 
varying m\n- This is shown explicitly by panels (c) and (d) of Fig. [31 which have the same 
mi/2 — 400 GeV as in panel (a) but with mo = 450, 550 GeV respectively. Going from panel 
(a) to panel (c), the decrease in m results in lower m^ and there is no longer any overlap 
between the region allowed by rrih and the region where the stop is the NLSP With higher 
mo as in panel (d), we get heavier ti and tit, and hence we need higher X t and hence Aq in 
order to approach t\ — x degeneracy. However, higher X t in turn has a problem with the 
rrih constraint as described above. This illustrates our conclusion that, indeed, a stop NLSP 
scenario is not possible within the CMSSM. 

6 Stop NLSP in the NUHM 

We next study the MSSM with non-universal Higgs masses (NUHM) [41]. In the NUHM, 
the soft supersymmetry-breaking masses in the Higgs sector, m\ and 7712, are not necessarily 
equal to the sfermion soft mass mo at the GUT scale. Using the radiative electroweak 
symmetry breaking conditions, we can characterize the new parameters as ix and itla (the 
pseudoscalar Higgs mass), both values being defined at the weak scale: 

„ m? + m 2 ,tan 2 /?+ ±(1 -tan 2 /?) + aL 1} , . 

fi = — 7^— — , 20 

tan 2 /3- l + A^ 

and 

m\{Q) = m 2 (Q)+m 2 (Q)+2/i 2 (Q) + A A (Q). (21) 

Many different sparticles could be the NLSP in this model, in different regions of the NUHM 
parameter space. These include the lightest neutralino x? the lighter stau fi, the selectron 
(smuon) e R (P-r), the lighter stop £1, the up squark (charm squark) u R (cr), and the tau 
sneutrino v T . Thanks to the RGE, the sbottom tends to be heavier than the stau, so unless 
we have non-universal soft scalar supersymmetry-breaking masses for squarks and sleptons 
at the GUT scale, sbottom could not be the NLSP. The up and charm squarks could be the 
NLSP for very large \\x\ and rriA [42], where the Higgs soft mass-squared m 2 or m\ become 
negative at the GUT scale. However, we discard this possibility, preferring to impose a GUT 
stability constraint [41]. We also do not consider the other possible NLSPs in this paper, 
and focus only on the stop. 

In order to obtain a light stop, small \p\ is preferred and, more weakly, large txia- In 
Fig. IU(a) we plot various masses as functions of \x for the fixed values tua = 1400 GeV, 
tan/5 = 10, A = 2100 GeV, mo = 500 GeV (in the region of the values studied previously 
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m 1/2 = 600, m = 500, m A = 1400, A = 2100, tan p = 10 m 1/2 = 600, m = 500, (I = 750, A = 2100, tan p = 10 
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Figure 4: Sparticle masses in the NUHM as functions of (a) \x with = 1400 GeV, and 
(b) m A with /I = 750 GeV, both with tan/? = 10, A = 2100 GeV, m = 500 GeV and 
m i/2 — 600 GeV. We plot the masses of the lightest neutralino (solid red), the second lightest 
neutralino (dashed red), the lighter stau (purple solid), and the lighter stop (orange solid). 
In the shaded region, the stop becomes tachyonic. The Higgs likelihood constraint is shown 
by the vertical dashed red line. We also draw a horizontal dashed line at 220 GeV to indicate 
the lower bound on the stop mass. 



within the CMSSM) and = 600 GeV (close to the minimum allowed so that m x > 
mf > 220 GeV). We see that there is a small range \x ~ 730 — 770 GeV where the stop 
is the NLSP and survives the phenomenological constraints. The stop relic density in this 
region is fl^ h 2 ~ 10 -4 , and the neutralino mass is about 250 GeV. We note that this region 
is nearly excluded by the Higgs mass constraint. Indeed, if we had taken the LEP limit on 
m h at its face value, this region would have been excluded. Panel (b) shows the spectrum 
as a function of m A , for fi = 750 GeV. Here we see that is a decreasing function of tua, 
and there is a region where a stop NLSP is allowed, at ~ 1350 — 1500 GeV. 

In Fig. Owe show spectra as functions of (a) m 1 / 2 , (b) m , (c) A Q and (d) tan/3. The 
lightest neutralino, which is bino-like in the cases shown, has a mass that depends essentially 
on mi/2 only. We see that increases as mm, ttlq or tan/? increases or Aq decreases, while 
rrif increases as m 1( / 2 , A Q or tan/? decreases or m increases. One might attempt to increase 
mi/2 to obtain a heavier NLSP. However, if we want a stop NLSP, in order to make the 
stop lighter than neutralino, we need to compensate the increase in mj by (say) increasing 
Aq, but the constraint on mi/2 would in turn render it more difficult to satisfy the mh 
constraint. We see allowed stop NLSP regions in panel (a) for m ly / 2 ~ 595 — 605 GeV, in 
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Figure 5: Sparticle masses in the NUHM as functions of (a) m 1( / 2 , (b) m , (c) A and (d) 
tan (3, with other parameter values specified in the legends. 
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m 1/2 = 600, m = 500, m A = 1600, A = 2100, tan p = 10 m 1/2 = 600, m = 500, |i = 800, A = 2100, tan p = 10 
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Figure 6: Same as Fig^ but with m t = 172.7 GeV. 



panel (b) for m ~ 490 - 510 GeV, in panel (c) for A ~ 2090 - 2110 GeV, and in panel (d) 
for tan (3 ~ 8 — 15. 

It is well known that the Higgs sector is sensitive to the top sector. Therefore, it is 
instructive compare with the result that would hold if mt = 172.7 GeV (an older experimental 
value of m t that is still within one a of the present central value). We also adjust so as 
to improve the overlap of the constraints, with the result shown in Fig. [6j We see that the 
Higgs likelihood constraint becomes as strong as the nominal rrih taken at face value, which 
is at n ~ 770 GeV, and that there is only a tiny region allowed. Thus, postulating a larger 
m t does not resolve the dilemma of the Higgs likelihood, even though the face value of is 
lifted. The allowed region of /i in panel (a) of Fig. [6] is narrower than in the corresponding 
panel of Fig. HJ and shifted to larger values of /i. Likewise in panel (b), comparison with the 
corresponding panel of Fig. H] shows that the allowed region is smaller and shifted to larger 
values of m^. 

We conclude that there are some small regions of the NUHM parameter space where the 
t\ is the NLSP, with a cosmological abundance that would correspond (in the absence of stop 
decays) to a relic density Qjh 2 ~ 10~ 4 . Typical allowed values of the NUHM parameters 
are m 1/2 ~ 600 GeV, m ~ 500 GeV, A ~ 2100 GeV, fj, ~ 750 GeV, m A ~ 1400 GeV and 
tan (3 ~ 10. We now discuss the cosmological evolution of such a scenario. 
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7 Cosmological Evolution of Metastable Stops 



We expect the metastable stop squarks and antisquarks density to have frozen out after 
coannihilation at a temperature Tp ~ ^ 7 GeV ^, when the age of the Universe t ~ 

10~ 9 s One might have expected a primordial stop-antistop asymmetry comparable to that 
for conventional baryons, but this would have been eradicated by stop-stop annihilations. 
The remnant stops and antistops would not have decayed before the next major event in 
standard Big Bang cosmology, namely the quark-hadron transition when t ~ 10 -6 to 10~ 5 s. 
At this point, they would have hadronized into sbaryons, antisbaryons and mesinos. 

Simulations and data from relativistic heavy-ion collisions indicate that the relative abun- 
dances of hadrons produced at the transition may be modelled by assuming an effective 
hadronic freeze-out temperature Tf ~ 170 MeV. We recall that all the heavier sbaryons 
would decay into the lightest state, which is charged, whereas the mesinos would all 
decay into the lightest T° state, which is relatively harmless. The mass difference between 
the A;, baryon and the B ±,Q mesons is about 365 MeV, and we expect a similar separation 
for the sbaryons and the T mesinos. Modulo spin and flavour counting factors, this 
would suggest a suppression by a factor ~ 10 for the abundance of the sbaryons and antis- 
baryons relative to the mesinos and antimesinos. This would suggest a priori an equivalent 
relic density Q\~h 2 ~ 10~ 5 , which might have been large enough to change significantly the 
subsequent abundances of light elements via bound-state effects. The S~ state, which is 
relatively long-lived, would have an abundance about a factor of 3 lower than that of A^. 

It has been argued recently that, following hadronization, stop hadrons hf would capture 
each other and form bound states, leading to the further annihilation of stops [27]. The rate 
for this process is approximately rih-(T '/ 'm^ x l 2 m~ 2 , where n^- — VLh~p c /mi and (T/mj) 1 / 2 
is the relative velocity in the capture process. The resulting abundance of stop hadrons is 
determined by equating this annihilation rate to the Hubble expansion rate: 



T 



n h ^ ( T\ 1/2 /8vriV T 2 



m 2 \mij V 3 Mf 

This yields a relic hf number density 



(22) 



n 1/2 

^^ a it P ^ < 23 > 



8 Note that the stop has a larger (co) annihilation cross section than does the neutralino, and hence 
would have smaller freeze out temperature relative to its mass compared to the neutralino with freeze-out 
temperature of m x /20 to m x /25. 

9 If the stop-neutralino mass difference is also ~ 5 GeV or smaller, which is quite possible in the allowed 
NUHM regions found in the previous Section, the stops would have been accompanied by a significant 
admixture of metastable neutralinos, which we discuss later. 
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where we use the estimate N = 0(50). The hf number density is minimized by annihilations 
at a temperature close to the formation temperature of the bound states, namely T ~ 200 
MeV, which yield 



T ~ 5 x 10~ 17 . (24) 



This corresponds to an effective flh f h 2 ~ 2 xl0~ 6 and = nf l ~mf l -/n 1 ~ 4 x 10~ 14 GeV 
for mi ~ 200 GeV. With this abundance, the late decays of stop hadrons would probably 
not cause problems with light-element abundances, even after allowing for uncertainties in 
the treatment of showers generated by hadronic decays 3 

However, the catalytic effects of bound states of the negatively-charged metastable relic 
antisbaryons A~ are potentially more dangerous, as emphasized above. Fortunately, as we 
now show, their abundance would have been suppressed following hadronization and before 
BBN, by annihilation with conventional baryons to produce T antimesinos and conventional 
mesons. The rate for the annihilation of heavy baryons is approximately n p m~ 2 . Since 
n p ^> and there is no velocity suppression for this annihilation process, since av — > 
constant as T — > 0, the rate for annihilation with nucleons would be much larger than the 
capture process [27] discussed in the previous paragraph. The abundance of antisbaryons 
would remain in equilibrium until the annihilation rate 



n-. 



1.5 x 10- 10 T 3 



* : , (25) 



where we use r\ = n v jn^ = 6.1 x 10 10 , became of the same order as the expansion rate 

^(OnPc) 172 1.2 x 10" 4 T 3 / 2 




M P Mf 



(26) 



where we use Q m ~ 0.27 for the matter density as a fraction of the closure density today. 
The rates (125]) and (126]) become comparable when 

/ 2 \ 2 /3 

r~10M^J ~0.1eV. (27) 

The A~ abundance would be suppressed by a factor exp(— 2 x 10 5 ) already by the time of 
the onset of BBN at T ~ 1 MeV, and the abundance of the dangerous A~ state would be 
driven to an extremely low value ~ (m^T) 3//2 e _m *' /T by the time of freeze-out (127]) of 
annihilations with relic baryons. 



10 We note in passing that the mechanism of [27] ceases to reduce the T abundance once the temperature 
falls below about 10 MeV, i.e., after about 1CP 3 s. 
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The relative abundance of the less dangerous A~ sbaryons would also have been strongly 
suppressed, as its annihilations with conventional antibaryons would have continued in equi- 
librium down to temperatures T ~ 20 MeV, when 



ml "V 3 M P ' 

where N — O(10) is the number of relativistic degrees of freedom at T ~ 20 MeV. 

Thus, during and after BBN, one must contend only with the neutral T mesinos and 
their antiparticles, with which they mix [30]. Their decays are relatively innocuous, as 
discussed above. In principle, they might also bind with deuterium to form superheavy 
nuclei and thus could change the Helium abundance. However, we would not expect the 
Helium abundance to be greatly affected: it is at the level of 24% of the baryon density and 
hence has Vtueh 2 ~ 0.01, whereas by the start of BBN the T mesino density would be down 
at the 10~ 6 level in the cases that we consider here, so the ratio of the number densities 
(assuming a ~ 200 GeV stop mass) would be (9(10 6 ). We note also that the T mesinos 
would react with conventional baryons to regenerate A~ sbaryons after their annihilations 
with antibaryons had frozen out, i.e., at T < 20 MeV. However, this process could generate 
a A~ density of at most 10~ 6 , which would not be problematic. 

This analysis therefore finds no cosmological problems with most of the small region of 
NUHM parameter space discussed in Section 6 where the stop is the NLSP. 



8 On the Possibility of a Metastable Neutralino 

The conclusion of the previous Section may, however, be modified if the stop NLSP is nearly 
degenerate with the lightest neutralino Xi a possibility suggested by our analysis of the 
NUHM parameter space. The \ ma y decay either into the stop (or antistop) and some other 
particles: x ~~ * ti + X, t\ + X', or into G + 7. The latter decay is very slow, being of 
gravitational strength [4]: 

^>=iii:|H)>^)' (29) 

where C X1 = (Oi x cos9 w + 2x sm9 w ). The decay rate depends on mg, and for mg ~ 
1 - 200 GeV we get T( X -> G7) ~ 10~ 28 - 10" 34 GeV corresponding to lifetime of ~ 
10 4 — 10 10 s. If the lightest neutralino is almost degenerate with the stop, the two-body decays 
X — > ti + 1, ti+t, which would be rapid but require a large mass difference m x — > m t , 
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are not available. As we discuss below, the decays x ~~ * ^i + A, t* + X' may be very 
suppressed if the x an d t\ are sufficiently degenerate. If they are suppressed sufficiently for 
the x lifetime to exceed about 10~ 5 s, neutralino relics would be present during and after 
hadronization, and the cosmology of the stops produced in their decays would differ from 
that of the thermally-produced stops discussed in the previous Section, as we show below. 

To show that the abundance of neutralino relics could indeed be significant, recall that 
when the stop and neutralino are almost degenerate, i.e., m x — vnq less than the freeze-out 
temperature ~ m^/30, neutralinos are kept in thermal equilibrium through co-annihilation 
processes [44] such as xX <-> t\X' where X, X' are Standard Model particles. We can 
approximate the neutralino relic density between the epochs of supersymmetric freeze-out 
and relic decays by 



tt x ~ - tt T ~ 0.25 tt NLSP , (30) 



1 

3 

modulo the Boltzmann factor from the mass difference. Here the factor 3 is due to colour, and 
^nlsp — ^ x + £\is the relic density one obtains from a standard coannihilation calculation, 
which leads for near-degenerate x an d h to flNisph 2 ~ 1.6 x 10~ 4 and hence fl x h 2 ~ 4 x 10~ 5 . 
Therefore the novel cosmology of the stops produced in neutralino decays is potentially 
important. 

If the x ~ ti rnass difference is larger than about 1 GeV, the neutralino may decay into 
three bodies, e.g., x ~~ > ti + s + 7r~, and into four bodies, e.g., x ~~ > ti + s + i + v and conjugate 
modes, and also into the corresponding final states with s — > d and/or i + v —>■ q + q. By 
analogy with the result of [43], the four-body semileptonic decay rate can be approximated 
as: 

T( X - t x Uv) = {0 G ^ Vtsl \ x [4|^ L | 2 (5m) 10 - 8Re[glg R }m t (SmY + 6|^| 2 m t 2 (5m) 8 )] 

(31) 

where we allow for both t\ and t* modes and a factor of 3 for colour. We have neglected the 
final-state fermion masses and use 

V 2 t L 3 cos 9w J rn w sm(3 J 

9r = ^^cos^t —O x4 -sm6 t - — O xl > , (33) 

V2 t m w smfj 3cosd w J 

where O is the neutralino diagonalization matrix T M^O = M^ ag . We see that the last 
term in the square bracket in (1311) is dominant, due to the large magnitude of m t . Therefore, 
we can approximate further to obtain 

T(x - hslv) ~ O(-) (2?r)5m , mf2 • (34) 
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For m^ = 220 GeV and 5m = 1 GeV, for example, and including the ev, \xv and du final 
states, we estimate T(x — > hsff) ~ (9(10~ 26 ) GeV, i.e., a x decay lifetime ~ 10 2 s. 

The partial decay time would be even longer for yet smaller 5m, but it would be necessary 
to take into account bound-state effects in the t^q and qq channels, and even in the full four- 
body final state. A naive scaling by (5m) 8 would yield a lifetime in excess of 10 10 s for 
5m ~ 0.1 GeV. It is therefore possible that the partial lifetime for x - ► h + X, t\ + X decays 
could exceed the partial lifetime for x ~^ G + 7 decays. 

As already remarked, if the x lifetime exceeds about 10~ 5 s, the t\ decay products ap- 
pear after the quark-hadron transition, and the discussion of stop cosmology given in the 
previous Section must be modified. If t x exceeds about 10" 3 s, the mechanism of [27] be- 
comes ineffective for reducing the abundance of stop hadrons produced in x decays, which 
therefore remain comparable to Q x h 2 ~ 4 x 10~ 5 , and also annihilations with antiprotons 
become ineffective for reducing the A~ abundance. On the other hand, annihilations with 
baryons remain effective for converting the dangerous A~ sbaryons into relatively innocuous 
T° mesinos for any x lifetime up to about 10 14 s. If r x < tj , the BBN/CMB constraints on 
electromagnetic and hadronic stop decays are insensitive to r x , and depend only on . If 
t x > r^, the the BBN/CMB constraints on electromagnetic and hadronic stop decays should 
be evaluated with tj replaced by r x . Finally, if the partial lifetime for x h + X, t\ + X 
decays exceeds that for x ~^ G + 7 decays, the x decay contribution to the relic stop density 
is diluted by the ratio of the partial decay rates. 

We therefore distinguish five metastable neutralino cases. 

1. If r x < 10~ 3 s, the residual suppression of the T° mesino density the mechanism of [27] 
may still be sufficient to evade the BBN/CMB constraints on electromagnetic and 
hadronic stop decays. 

2. If 10~ 3 s < t x < 10 4 s, the the mechanism of [27] is ineffective, and the indirectly- 
produced T° mesinos still have Vl x h 2 ~ 4 x 10~ 5 . Since the neutralinos decay before the 
epoch when the light-element abundances constrain the electromagnetic and hadronic 
stop decays, these depend on the value of r$ in the usual way. The relic density of T° 
mesinos respects the limit obtained in [16] by considering electromagnetic decays, but 
not the hadronic limit given in [15], which will be re-evaluated in [45]. 

3. The usual light-element abundance constraints also apply if 10 4 s < r x < r-jr , and the 
relic density of T° mesinos is again marginal. 

4. If t x > ti and 10 4 s simultaneously, the usual light-element abundance constraints 
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still apply, but they should be evaluated with tj- replaced by r x + . In addition, we 
should consider also the effect of neutralino decay itself to the light element abundance. 

5. In both the last two cases, the scenario may survive more easily if the partial lifetime 
for x ~^ ti + X, t* + X decays exceeds that for \ ~^ G + 7 decays, in which case the 
abundance of T° mesinos is suppressed by the ratio of the partial decay rates. 

9 Summary 

The scenario of gravitino dark matter with a stop NLSP has rich phenomenology, for both col- 
lider experiments and cosmology. Unfortunately this scenario is disfavoured in the CMSSM 
and NUHM, due to the existing collider limits and the close relation between the masses 
of the light stop and the light Higgs boson. However, we find that there could still be a 
small allowed region, at least within the NUHM, with wiiw ~ 600 GeV, mo ~ 500 GeV, 
A ~ 2100 GeV, n ~ 750 GeV, m A ~ 1400 GeV and tan/3 ~ 10. 

Much of the discussion in this paper on the stop NLSP scenario could be applied as well 
to a scharm or sup NLSP in more general MSSM models with a gravitino LSP. In such a case, 
we expect that the Higgs constraint would be less stringent. Presumably the late production 
of the NLSP would be less important since, with the light u, c quarks in place of the heavy 
t quark, two-body decay channels would probably be available, unless the mass difference 
between neutralino and the NLSP is very small (less than m UjC ), in which case the neutralino 
would decay directly into the gravitino. The case of a sbottom NLSP would in this respect 
be intermediate between the stop and scharm cases. However, as noted above, the lightest 
mesino would, in this case, probably be charged, and hence capable of bound-state catalysis 
of dangerous light-element transmutations. 

For a possible realization of models with a sup NLSP, see [42], where one could get light 
scharm and sup in the NUHM by assuming very high values for /i and tjia, which would 
mean violating the GUT stability constraint. Another possiblity would be to postulate 
non-universality between the first two generations and the third, and between squarks and 
sleptons. It is interesting to note that the suppression of the lighter squark eigenmass via 
diagonalization is much less for sup than for stop, so it is not necessary to postulate large 
Ao in these scenarios. 

As Hamlet said: 'There are more things in Heaven and Earth, Horatio, than are dreamt 
of in your philosophy' This comment certainly applies to supersymmetric phenomenology. 
There are surely many important aspects of the stop NLSP scenario that we have overlooked 
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in this paper, and many other NLSP candidates could be envisaged, beyond the stop, the 
neutralino and the stau, which are the three options usually considered. If supersymmetry 
is the 'surprise' most expected at the LHC, one should perhaps expect it to appear in an 
unexpected way. The dominant signature might not be the 'expected' missing transverse 
energy, but rather some brand of metastable charged particle, which could well have strong 
interactions, like the stop considered here. 
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Appendix: Three-Body Stop Decays 

We present in this Appendix some details of the calculation of the three-body stop decay, 
using the following notations. For the top contribution, we introduce: 



A t = o 



B; EE - 



(cos 6 1 + sin%), 
[cos 9 1 — sin 9i). 



(35) 



(36) 



and for the sbottom contribution we introduce: 



Oj ee (sin^,cos%), 
hi ee (cos 9~ b , -sin 9~ b ). 



(37) 
(38) 



We also define 



Ki EE (COS ^ COS ^, — COS COS ^). 



(39) 



In the chargino contribution we have: 



Vi ee -(V i2 sm/3 + U i2 cos(3), 
A, t ee hv i2 smP - U i2 cosP), 



(40) 



(41) 
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and for the low-to-moderate range of tan (3 we have: 

oc i I 92-mt sin L sin 6~ t 

2Si = -# 2 cos0 L + -= — — (42) 

V Zmw sin p 

op j, sin0 L sin fl f - 

2Pi = -# 2 cos0 L 7 = — — (43) 

\J lm w sin p 

where cos0L,±sin0£ are elements of the matrix V that diagonalizes the chargino mass 
matrix, and expressions for S2 and P 2 may be obtained by replacing cos0l — > — sin0L and 
sin0£, — > cos0l in the last equations. 

For the square of the top contribution to the decay amplitude, we have: 

(44) 



where 

/ii = m\m\-fl, (45) 

^2 = -2/ 1 /3 + 2/ 3 2 + m 2 w (-/ 2 + 3/3 + m|), (46) 

fea = -2m 2 4 ,g 1 2 -2(/ 1 -/ 3 )(2/i- 2/ 3 - 3m^). (47) 

From the energies Pi = p? and E w = k°, we can define variables x,y by P x = m~ x/2 and 
E w = rri~jj/2. In turn, this allows to express all the inner products of momenta that appear 
in the functions Wu in terms of x and y. Then, q\ = m-^l + r\ — x) and the functions fa 
are given as follows: 

/1 = -fy, (48) 

2 

771- 

h = -fx, (49) 

2 

Tfl~ 

h = -f(-l-n-r 2 + x + y). (50) 

where r 1 = m 2 ~/m 2 ^ r 2 = m^/m^. 

Similar expressions can be written for the square of the sbottom contribution to the decay 
amplitude, namely: 

= (51) 



where: 
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and: q\ = m| (1 + r 2 - y), q 2 .pi = /2 - h and p x .p 2 = f 2 - h~ m 2 ~. 
The square of the chargino contribution takes the form: 



l^x+l 2 = E^ P x>fe)^ + fe)^ x ; (53) 

where 

32 

W X + X + = o^2 ^2 ( 2m W m2 G + /I) [ 2 /4( S (l)ij^G m X + S (2)ii/s) + Pl-P2 (™ x S (3)u ~ ?3 S (2)ij)] 

(54) 

and ql = m| + + 2/ 3 , / 4 = (m|/2)(2 -x -y), f 5 = / 3 + m|. Also, we define = 

(AiAj-ViVj^SiSj+PiPj), z m = (AA+^)(^+^0-^^^ 

and E (3)ij = (AA + + p i p j) + (A^ + ViA^SiPj + P^). 

The interference between the top and chargino contributions leads to the following ex- 
pression for W tx +: 



2 2m + 

W txf = 2(T H - m t T 2t ) + 2m xf (T 3t - m t T 4i ) - - (T 5l - m t T 6l ) - (T 7t - m t T 8t ) (55) 



The functions Tj may be written as follows: 

T « = ^^[-2f!ml + (f 2 fi-my i f 3 )h 4 -2f 1 f 2 f 3 (-f 2 + 2f 3 + ml) 



-2flm 2 w + f 2 m 2 w (fi + 2m 2 d ml) - hm 2 w {2^ + f 2 f 3 + m z ~f 3 - 2m%m\ 
+m 2 w (f 2 2 - m 2 d ml)h 5 + tf(2f 2 f 3 - 2m|/ 2 + m|(4/ 3 + 2m| + m 2 w ))], (56) 



T 2 , = 4^[(/i-/3)(2/3 + m|)(m|/ 1 -/ 2 / 3 ) 

+m 2 w (-f 2 f 2 3 - /|(/ s + m|) + m%{hh + h™\ + m|m?))], (57) 
T 3l = -^[2A 2 m|(-|-/ 2 ) + ^3 2 H-<) 

+/i/ 3 (2/ 2 2 - 2/ 2 m| - m% + m|m|) 

-m^/f - m|m?)(2/ 2 - / 3 - 2m|)], (58) 



25 



r : - -^^[(/i-/3)(^|/i-/2/3)+m 2 w (/|-m|m|)], (59) 



T 5l = ^^[-2/ 1 2 (/ 2 -m|)(/3 + m|) + /3 2 / 2 (2/3 + m|)-/ 3 2 m? i (/2 + 2/3) 

m G m W 

+hh(2f 2 - 2m|/ 2 - (2/ 3 + m|)(m| - m?)) 
+m^/ 2 (-2/ 2 + 5/ 2 / 3 + / 3 2 + 2m|/ 2 ) 
-m 2 y m|(2/ 2 (/ 3 - m|) + (/ 3 + m|)(/ 3 + 2m|)) 

mUMfs ~ 2m|) - m|(/ 3 - 2m|)) + 2m 4 w (/ 2 - m|m? )], (60) 

' n Q ln W 

+/i/3(2/ 2 / 3 + ^|(3/ 2 + 2/3) + m|) 

+h{-m 2 w fj - m\m\{h + - / 2 / 3 (2/ 3 + m| + m 2 w ))], (61) 
^ = 4^[A 2 (-2^|/ 2 + 2mi) + / 2 / 2 (m|-m 2 i ) 

+/i/ 3 (2/| - 2/ 2 m| - m| + m|m|), 

-/ 2 m 2 y (/ 2 (2/ 2 - / 3 ) - 2m|(/ 2 + / 3 )) - m\m\m 2 w {-2j 2 + 3/ 3 + 2m|)], (62) 
T 8l = ^^[2m!/ 2 + /3(/ 2 /3 + /3<+2/^ 

m G m W 

where a u = A 3i A~ t + A 4i B h a 2i = A 3i Ai - A 4i B h a 3i = E 3i A t ~ + E 4i B h a 4i = - T, 4i B h 

with A 3i = - AiPi, S 3i = VJSj + APi and A 4l = - AS h S 4l = V.P, + 4$. We 
also denote h 4 = 2/ 3 + m% — m~ , h 5 = 3/ 3 + 2m~ + mL. 

The interference between the top and sbottom contributions, W t ~ b may be written as: 

^ ^^2 ^/l ^^2 ^ ^ 

^ = 2 2 " * 2 2' 2 G ["4*1* - A^ifci] + [m^rapi - /ii2w] , (64) 
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where 



Y lpi = AZ ll [{h + f 2 )m 2 d + {f 2 -h-2m 2 d )m 



G 1 \ JZ Ja G> t\\ 

4 Z 2i [(fi + f 2 - mlJmtmQ], (65) 
Y lkt = 4Z H [2/ 1 (/ 2 -m|) + / 3 (m|-m? i )+m^(-/ 2 + m|)] 

-4 + / 3 - m 2 w )m t m~\, (66) 

= -^[/3 3 (^| + + /2m|(/ 3 2 + 2/ 3 m|-m|(m| + m?)) 

+m|m|(-/ 3 2 + / 3 m| + m|(2m| + m^)) 

-fliH^l + 2m V + m U 2m i + m ^)) 

+m|/i(/ 2 (/ 2 + / 3 + m|) + m|(/ 2 - / 3 - 3m|))] 
+/ 2 2 (/s + m? + m^) + m|(/ 3 2 - m|(m? + m^)) 

+m d m t f 1 (-f 2 (f 2 + / 3 + m|) + 2m|m|)], (67) 

*2ki = -^[(/|-/2/3-m|m t ? 1 )(/3(-m| + m|)+m 2 4 ,(/ 2 -m|)) 



-/!(2/ 2 3 - 2/ 2 (/ 3 + m|) + m|/ 3 (m| - m|) 

+m|/ 2 (2/ 3 - 2m| - m^) + m|(2m| + m^)) 

-2m|/ 2 (/ 2 -m|)] 
4 Z 

+ — ^[-hmQm t {-fl + (/ 3 + m~)m| + / 2 (/ 3 + m 2 w )) 
m~ 

+fim d m t (-f 2 + 3/ 2 / 3 + m|(/ 3 + m| + m 2 w ) - 2.Am|)], (68) 



and Z u = AfL 5i + BfA 5i , Z 2i = A f T, 5i - BtA 5i , S 5i = (a i + h i )/2, and A 5i = (a; - 6j)/2. 
Finally, the interference between the chargino and sbottom contributions is given by: 

2 

= 2 [ Xi p^ ~~ m x + X 2 P ij] - - [JCiwj - m x +X 2kij \ (69) 
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where 

Xirt = ~^-[(-A-A + 2/3 + ^1 + ^)^6 + h 7 -myi)\, (70) 

m G m W 

Xmi = -^^[A 2 (/3(/2 + /3)+m|(/ 2 -2/ 3 )-ml) 

-fihUl - / 2 (3/ 3 + m%- 2m 2 w ) + / 3 (2/ 3 + m| + m| + m^)) 
+™wUl - / 2 2 (3/ 3 + m|) + (/ 3 + m|) 2 m| + / 2 (2/ 3 2 - m|m| + / 3 (m| + m^)))], 

(71) 

X 2TO - = -4^V(-/ 2 + /3 + m|)[/ i6 + / i7 -/ 1 2 m|], (72) 
X 2fc{J = -^V(-/ 2 + /3 + m|)[/i 6 + / i7 -/ 1 2 m|], (73) 

and /ib = /i/ 3 (/ 2 - / 3 ), ^7 = m^(/ 2 (-/ 2 + /a) + ™~™|), and cr Hi = r ay T^ - T 2ij A i , a 2ij = 
TujVi + T 2 ijAi. We also define r\, = E 5j Si + A 5j Pi, T 2i j = S 5 j-Pi + AsjSi, with E 5j , A 5j 
defined previously. 
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